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Ultrasonic  measurements  of  wave- propagation  characteristics  in  composite  and  heterogeneous 
materials  provide  an  excellent  means  to  study  their  mechanical  properties.  In  recent  years  we  have 
studied,  both  theoretically  and  experimentally,  characteristics  of  elastic-wave  propagation  in  particle- 
reinforced  composites  and  heterogeneous  materials  as  well  as  in  homogeneous  and  laminated  fiber- 
reinforced  composites.  Comparison  of  theoretical  predictions  with  obervations  of  wave  velocities 
has  shown  good  agreement  and  has  provided  a  way  to  evaluate  microstructural  dependence  of 
mechanical  properties  of  these  materials.  Modeling  predictions  coupled  with  observations  can  also 
be  used  to  obtain  mechanical  properties  of  the  reinforcing  phase,  which  are  sometimes  not  easily 
obtained.  In  this  paper  we  present  results  of  some  of  these  recent  studies. 


We  also  present  results  of  our  study  of  changes  in  phase  velocities  and  attenuation  caused 
by  interface  layers  between  the  reinforcing  phase  and  the  matrix.  We  show  that  this  third  phase 
measurably  modifies  the  dispersion  behavior.  This  should  lead  to  effective  characterization  of 
interface  layer  properties  by  ultrasonic  methods.  ^ 

INTRODUCTION 


Determination  of  effective  elastic  moduli  and  damping  properties  of  a  heterogeneous  or 
composite  material  by  using  elastic  waves  (propagating  and  standing)  is  very  effective.  Several 
theoretical  studies,■,,  show  that  for  long  wavelengths  one  can  calculate  the  effective  wave  speeds  of 
plane-longitudinal  and  plane-shear  waves  through  such  a  material.  At  long  wavelengths,  wave 
speeds  thus  calculated  are  nondispersive;  they  provide  values  for  the  static  effective  elastic 
properties.  References  to  other  studies  can  be  found  in  those  cited. 

We  present  results  of  some  of  our  recent  studies  of  phase  velocity  and  attenuation  of 
plane  longitudinal  and  shear  waves  propagating  in  a  medium  with  microstructurc.  Microstructures 
studied  were  either  inclusions  or  fibers.  In  the  first  case,  we  examined  the  effect  of  inclusion 
shape,  volume  fraction,  and  elastic  properties  on  wave  speeds.  We  studied  inclusions  with  their 
interface  layers  separating  them  from  the  matrix  medium.  For  fiber-reinforced  materials,  we 
studied  continuous  aligned  fibers.  In  this  case,  the  medium  behaves  anisotropically  because  of  the 
alignment  of  the  fibers. 


The  theoretical  model  for  these  microstructural  studies  used  a  wave-scattering  approach  and 


predicted  the  macroscopic  isotropic  elastic  properties  for  the  case  of  random  orientation  of 
inclusions  and  for  spherical  inclusions.  For  aligned  fiber-reinforced  materials  the  model  gives  the 
anisotropic  elastic  properties.  The  scattering  approach  led  also  to  an  estimation  of  attenuation  and 
dispersion  of  waves  via  the  optical  theorem. 


The  scattering  approach  applies  to  a  macroscopicallv  homogeneous  medium,  infinite  in 
extent.  For  a  bounded  medium  such  as  laminated-plate  structures  we,  developed  a  hybrid 

numerical  technique  to  analyze  dispersion  of  guided  waves.  In  this  paper  we  present  some  of  our 
computational  results  showing  the  effect  on  dispersement  characteristics  caused  by  interface  soft 
layers  between  stiff  layers  in  a  p'ate. 


The  experimental  methods'0-1*-1*  consisted  of  a  pulse-echo  technique  and  a  resonance- 
method  These  were  chosen  to  provide  the  advantages  of  small  specimens  and  low  inaccuracy 
For  details  of  the  experimental  techniques,  the  reader  is  referred  to  the  references  cited 

THEORY 


Scattering  by  a  Single  Inclusion 

Consider  a  single  elastic  ellipsoidal  inclusion  with  material  properties  X’,  it'.  p‘  embedded 
in  an  clastic  matrix  with  properties  X.  it.  p.  Assume  that  the  inclusion  is  separated  from  the 
matrix  by  a  thin  layer  of  elastic  material  with  properties  X,.  itx.  P,.  which  are  variable  through  the 
thickness.  Here  X.  ft  arc  the  Lame  constants  and  p  is  the  density.  The  geometry  of  he  problem 
is  shown  in  Figure  l 

We  need  to  find  the  field  scattered  by  this  inclusion  when  it  is  excited  by  an  incident 
clastic  wave  A  low-frequency  approximate  solution  to  this  problem  was  presented  earlier'6  for 
the  case  of  no  intermediate  layer  between  the  inclusion  and  the  matrix.  It  can  be  shown  that  if 
the  exciting  field  is  given  by.  dropping  the  time  factor  c,<Jt . 
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Here  L  N  ***,  and  arc  spherical  vector  wave  functions'7  that  are  finite  at  r  -  0  and 
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where  v0  is  the  volume  of  the  ellipsoid  and 
e  =  k,C.  k,  =  cj/C,.  t  =  C,/C2 
3rJ,  k=0,2 
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C,  and  Cj  are  the  longitudinal  and  shear  wave  speeds  in  the  matrix  medium.  Expressions  for 

*  "••'i 

T  .  can  be  found  elsewhere’ 

MM 

Effective  Properties  of  a  Composite  Medium  with  Inclusions 

Once  the  scattered  field  caused  by  a  single  inclusion  is  known,  multiple  scattering  from  a 
number  of  inclusions  can  be  easily  calculated.  In  particular,  for  a  random  homogeneous 
distribution  of  ellipsosidal  inclusions 
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Here  c  =  - - .  n0  being  the  number  density  of  inclusions.  The  constants  P0  and  1%  depend 

3irabcn0 

on  the  geometrical  properties  of  the  inclusions  as  well  as  on  the  mechanical  properties  (X.m.X'.m') 
of  the  matrix  and  the  inclusions.  P,  is  simply  given  by 

i»  .  p'/p-I 


In  deriving  (4)  and  (5)  it  has  been  assumed  that  the  inclusions  are  similar  in  shape,  size  and 
phvMcal  properties. 

When  there  is  an  intermediate  layer  between  the  inclusion  and  the  matrix,  then  the 
coefficients  A^„  and  13MI,  cannot  be  obtained  in  exact  form;  they  have  to  be  determined 
numerically  for  general  elliposoidal  inclusions.  For  spherical  inclusions,  on  the  other  hand,  exact 
expressions  for  and  I3„„.  valid  at  arbitrary  frequencies,  can  be  obtained  if  the  intermediate 

layer  has  constant  properties.  Also,  for  spherical  inclusions  with  thin  interface  layers  having 
variable  properties.  A„„  and  can  be  obtained  approximately  if  h/X<<l.  h/o <  <  1 .  Here  h  is 

the  thickness  of  the  interface  layer,  a  the  radius  of  the  spherical  inclusions,  and  X  the  wavelength 
of  the  incident  wave.  Even  then,  the  problem  of  determining  phase  velocities  of  longitudinal  and 
shear  waves  in  the  composite  material  at  finite  frequencies  is  complicated  .  The  problem  is 
simpler  if  the  volume  concentration  of  inclusions  is  small.  This  is  discussed  in  the  following 

Consider  a  spherical  inclusion  of  radius  a  with  an  interface  layer  (Figure  2).  let  the 
incident  field  be  given  by 

i  ik.z  ik,z  ,,, 

u‘=e  e2+e  ex  (6) 


/here  k  =  w/C,.  The  scattered  field  is  then  given  by 
n  ^  =  u  P  +  us 


where  superscripts  p  and  s  refer  to  longitudinal-wave  and  shear-wave  components,  respectively. 
Now  it  will  be  assumed  that  within  the  interface  layer  the  elastic  coefficients  X,  and  m,  vary  with 
the  distance  from  the  center  of  the  inclusions  as 

X,(r)  +  2/x,(r)  =*  (X,'  +  2m,’)  f (r) »  a<r<a+h,  (8) 

M,(r)  =  m,'  g(r) ,  a<r<a  +  h, 

where  f(r)  and  g(r)  are  general  integrable  functions  of  r.  If  it  is  further  assumed  that  h/a <  <  1 . 
h/X<<l,  then  it  can  be  shown  that  the  displacement  components  satisfy  the  approximate  boundary 
conditions  on  r  -  a. 
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Superscript  (t)  refers  to  the  field  quantities  within  the  inclusion.  is  the  stress  tensor. 

To  this  order  of  approximation  the  traction  components  rrr.  rrg  and  are  continuous  at  r  = 
a  These  simplified  boundary  conditions  allow  the  single  scattering  problem  to  be  solved 
exactly 

When  r  -  oo.  one  obiains  from  Eq.  (7)  the  far-field  behavior  of  u S  when  the  incident 

wave  is  given  by  a  plane-longitudinal  wave  (the  first  term  on  the  right-hand  side  of  Eq.  (6)). 
Then  it  can  be  shown  that 
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For  the  incident  plane-shear  wave  (the  second  term  on  the  right  hand  side  of  Eq.  (6)1  one  finds 


u  P  ~  gP  (5.0)  e  r 


ikjr  ik2r 

u  s  ~  hs ,(0.0)  - -  c  g  +  h,s(0.$)  - -  e 


The  expressions  for  the  amplitude  functions  gP.  hP.  and  so  on,  can  be  obtained  from 
our  earlier  study".  Using  Eqs.  (10)  and  (11)  and  the  forward-scattering  theorem"  we  then  obtain 
the  equations  for  the  effective  wave  numbers  in  the  composite  medium: 

k  *a 

~-~r  -  1  +  — y  n0  gP(0)  (12) 

k,  k, 

for  the  longitudinal  wave  and 
k  *J 

— —p  *  1  ♦  h,s(0,0)  (13) 

k,  k, 

for  the  shear  wave  Since  gP(0)  and  h,s (0.0)  are  complex,  effective  waves  will  be  both  dispersive 
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arul  alternative. 


Effective  Properties  of  a  Fiber- Reinforced  Composite 

The  analysis  presented  above  for  inclusions  can  be  applied  also  to  a  medium  reinforced  by 

aligned  continuous  fibers.  In  that  case,  one  can  derive  equations  similar  to  (4)  and  (5)  for 

longitudinal  and  shear  waves  propagating  perpendicular  to  the  fibers.  Taking  the  x,  axis  along 

20 

the  fibers,  and  assuming  them  to  be  transversely  isotropic  about  this  axis,  it  was  shown  that 
for  Sll  waves  polarized  along  the  fibers  the  effective  wave  number.  S'.  is 

0  '  _  p  '  ( l-c(m- 1)  (m~  1 ) )  f. 

|T  p  (l+c(m-l)/(m+l)) 

where  p *  is  the  effective  density  and  m  =  C,,/p.  For  longitudinal  and  shear  waves  polarized  in 
a  plane  perpendicular  to  the  fibers  we  get 


p*(UcP0)[UcP2(UrJ)l 
1  -cPj(  1  -  T)  -  2c:P3Pj 
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Note  that  Cu.  C„,  _CJ4,  CM.  and  C,,  are  the  five  independent  elastic  constants  characterizing  the 
fibers  and  p  =  p[Uc(p'/p-l)]-  P0  and  Pa  are  defined  by  the  relationships 

K't  -  (X*m)  p(C«-m) 

P  =  - - i -  P  =  -  - - -  (17) 

0  K’t+m  1  J  C  ' 


where  K'-p  is  the  plane-strain  bulk  modulus  of  the  fibers.  The  remaining  two  elastic  moduli  of 
the  composite  are  obtained  from  relationships  derived  by  Hill  . 

Dispersion  of  Guided  Waves  in  a  Laminated  Plate  with  Interface  Layers 

In  composite  materials,  interfaces  between  the  different  constitutents  play  an  important  role 
in  determining  their  mechanical  behavior.  As  discussed  above,  the  effect  of  interfacial  layers 
between  inclusions  and  matrix  medium  on  composite  properties  can  be  evaluated  approximately 
using  a  scattering  approach.  For  laminated -composite  plate  structures,  the  effect  on  the  dispersion 
characteristics  caused  by  interfacial  layers  between  the  laminae  can  be  analyzed  in  detail.  The 
changes  in  dispersion  characteristics  should  provide  a  way  to  ultrasonically  evaluate  the  interface 
properties 

In  this  section  we  briefly  outline  a  theoretical  technique  to  study  guided-wave  propagation 
in  a  laminated  plate  with  interface  layers  For  simplicity  of  analysis,  we  consider  only  isotropic 
laminae  and  isotropic  interface  layers  Equations  governing  guided-wave  propagation  in  such  a 


plate  can  be  solved  exactly.  However,  if  the  laminae  are  isotropic,  then  an  exact  analysis  is 
extremely  complicated.  To  avoid  the  difficulties  associated  with  an  exact  analysis,  we  developed  a 
stiffness  method  in  which  each  lamina  (and  interface  layer)  is  divided  into  several  sublayers 
Polynomial  interpolation  functions  for  through-th.ckness  variation  of  displacements  in  each  layer 
are  assumed  The  interpolation  functions  involve  a  discrete  number  of  generalized  coordinates  that 

arc  functions  of  the  in-plane  coordinates  x  and  y  (Figure  3)  and  time  t  The  generalized 

coordinates  are  the  displacements  and  tractions  at  the  interlaces  between  the  adjoining  sublayers, 
thus  ensuring  continuity  of  these  quantities  at  the  interfaces.  13y  applying  Hamilton's  principle  the 
dispersion  equation  is  obtained  as  a  standard  algebraic  eigenvalue  problem  whose  solutions  yield 
the  dispersion  relations  and  the  \anation  of  stresses  and  displacements  through  the  plate  thickness 
of  the  plate.  This  technique  was  used  earlier”  to  study  wave  propagation  in  an  infinite 

periodically  laminated  medium.  The  application  of  the  technique  to  a  sandwich  plate  was  discussed 
in  detail  elsewhere”.  Here,  for  brevity,  we  present  only  the  results  for  a  particular  sandwich 
plate. 

RESULTS  AND  DISCUSSION 

Cast-Iron  Elastic  Constants 

Graphite-particle  shape  affects  cast-iron’s  properties,  both  physical  and  mechanical.  Several 
studies24-”  dealt  with  the  various  properties  of  cast-iron  as  they  depend  on  graphite  particle 

shape.  But  all  these  studies,  experimental  and  theoretical,  dealt  only  with  limiting  shapes:  sphere, 
rod,  and  disc.  They  failed  to  deal  with  arbitrary  aspect  ratio,  c/a.  of  the  particles.  Using  Eqs. 

(4)  and  (5),  we  calculated  Young's  modulus.  E,  of  cast-iron  for  various  value*  of  c/a.  The 
results  are  shown  in  Figure  4  along  with  the  experimental  observations  reported  by  various 

investigators  21-26.  In  this  figure,  calculated  results  are  for  two  different  volume  fractions-- 10  and 
12  pcrcent--a  range  that  contains  most  of  the  studied  cast-irons. 

The  two  upper,  nearly  horizontal,  curves  correspond  to  graphite's  upper  third-order 

elastic-constant  bounds.  The  two  lower  curves  correspond  to  the  lower  bounds.  From 

monocrystal  elastic  constants,  using  equations  by  Kroner  and  Koch28,  for  graphite,  Wawra.  et  al2’ 
calculated  third-order  elastic-constant  bounds.  They  found  the  following  effective  quasi-isotropic 
elastic  constants  E'=  4  17(1.34)  GPa;  m'=  1.41(0.45)  GPa.  Values  outside  parentheses  denote 
upper  bounds;  those  inside  denote  lower.  For  the  matrix  phase  we  took  the  constants  for  alpha 
iron:  E  =  206  GPa,  n  =  80.0  GPa. 

Corresponding  to  observation,  our  model  predicts  a  strong  dependence  of  Young's 

modulus  on  aspect  ratio.  Near  the  spherical  limit  (c/a  =  1).  E*  .  varies  slowly  with  c/a.  Near 
the  oblate-disc  limit  (c/a  *  0),  E*  varies  rapidly  with  c/a.  An  interesting  result  is  that  graphite's 
lower-bound  quasi-isotropic  elastic  constants  fit  observation  so  well 

Elastic  Constants  of  Graphite- Aluminum  Composite 

Figure  5  shows  the  microstructure  of  graphile-fiber-reinforced  aluminum  obtained  from  a 
commerical  source.  By  Archimedes’  method,  we  found  the  mass  density  of  the  composite  to  be 
2.013  g/cm*.  For  a  fiber  volume  fraction  of  0.70,  using  2.6523  for  aluminum,  the  graphite  fiber 
density  is  predicted  to  be  1.738,  very  close  to  the  manufacturer's  estimate  of  1.76. 

We  determined  the  nine  Cjj  by  measuring  eighteen  sound  speeds  on  four  specimen 
geometries  described  previously”.  For  brevity,  we  omit  further  description,  except  for  a  few 
salient  details:  bond--phenylsalicylate;  transducers--quartz.  x-cut  and  ac-cut;  frequcncies--5  to  6 
MHz;  specimen  size--16-mm  cube,  or  smaller  depending  on  specimen  geometry. 
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Tabic  I  shows  the  study's  principal  results.  Column  1  lists  various  elastic  constants. 
Column  2  gives  a  set  of  fiber  elastic  constants"-’2  We  chose  these  because  E„  agrees  closely 
with  the  E„  for  the  present  fiber.  Column  3  shows  measured  results.  From  the  measured 
results  and  the  predictions  from  Eqs.  ( 14) - ( 16)  together  with  Hill's  relationships,  we  predicted  the 
fiber  properties  shown  in  column  4  We  used  the  calculational  sequence:  C44,  Cbtl.  C,,.-CM.  »«„, 
and  Eu.  Column  5  shows  the  calculated  C,j.  Finally,  column  6  shows  the  ratio  of  column  5  to 
column  3 

Results  of  column  3  of  Table  I  show  that  the  studied  composite  shows  orthotropic  elastic 
symmetry,  which  is  approximately  transversely  isotropic  with  the  x,-axis  as  the  symmetry  axis. 
The  nucrostructure  in  Figure  5  also  suggests  transverse-isotropic  symmetry. 

Concerning  the  first-guess  graphite-fiber  elastic-constant  calculations,  we  find  reasonable 
agreement  for  C,,.  C,2.  C„.  and  C6S.  Thus,  the  criterion  of  choosing  a  graphite  elastic-constant 
set  based  on  En.  the  axial  Young's  modulus,  succeeds  partially. 


<■ 


One  can  obtain  a  better,  complete  graphite  elastic-constant  set  by  using  the  mode! 
equations  inversely  This  is  shown  in  column  4  and  it  is  seen  that  the  constants  appearing  in 
this  column  differ  significantly  from  the  first-guess  values  in  column  2. 

For  the  graphite-aluminum  composite.  Figure  6  shows  the  variation  with  temperature  of 
the  principal  elastic  stiffnesses  C y  (i  =  1.6).  For  comparison,  the  figure  also  shows  the 
temperature  variation  of  the  longitudinal  and  shear  moduli  of  the  aluminum  matrix.  The 
graphite-aluminum  Cjj-versus-temperature  results  in  Figure  6  show  strong  anisotropy.  The  largest 
change  occurs  in  Cu.  the  longitudinal  elastic  stiffness  perpendicular  to  the  fibers.  The  smallest 
change  occurs  in  C„.  the  longitudinal  elastic  stiffness  parallel  to  the  fibers.  These  changes  agree 
with  the  well-known  high  axial  elastic  stiffness  and  low  axial  thermal  expansivity  of  graphite 
fibers.  The  three  shear  moduli--C44,  C„.  C51--fall  between  these  extremes.  Among  the  shear 
moduli.  shows  the  largest  change;  this  reflects  the  low  C66  values  for  both  the  fiber  and  the 
matrix.  Almost  universally,  lower  elastic  stiffness.  C,  means  higher  dC/dT.  In  Figure  6,  the 
near  equivalence  of  dC44/dT  and  dC„/dT  reflects  the  approximate  transverse  isotropy  of  this 
composite. 

Interface  Effects  on  Damping  and  Phase  Velocities  in  an  SiC-ParticIc  Reinforced  Aluminum 

Equations  (12)  and  (13)  provide  implicit  relations  for  the  complex  wave  numbers  k,*  and 
kj*  if  it  is  assumed  that  the  inclusion  is  placed  in  a  composite  medium  with  the  effective 

unknown  dynamic  properties.  Then  k,2  and  k,2,  appearing  on  the  right-hand  sides  of  Eqs.  (12) 

and  (13)  will  be  replaced  by  k,2  and  k22  respectively.  These  equations  then  can  be  solved 

iteratively  for  the  unknown  k,*  and  k2* .  These  results  are  shown  in  Figures  7-10.  Note  that 
Im(k*)  measures  damping  and  Relk^)  measures  the  ratio  of  the  phase  velocities  in  the  matrix 
and  the  composite.  In  these  calculations  the  interface  layer  properties  were  assumed  to  vary 
linearly  across  the  layer  thickness  from  the  properties  of  the  particles  to  those  of  the  matrix. 

The  elastic  properties  of  the  particles  and  the  matrix  were  taken  as:  X'+2 m’  ■  4.742x10"  N/m  . 
It'  =  1. 881x10"  N/m’.  p'  =  3.181  g/cni’.  X+2 n  «  1. 105x  10"  N/m’.  n  *  0.267x10"  N/m*.  f  -  2.705 
g/cm  .  Calculations  were  performed  with  or  without  interface  layers  and  at  two  volume  fractions: 
c  ■  0.05.  0.15.  The  presence  of  the  interface  layer  decreases  the  damping  as  well  as  the  phase 
velocities. 

Dispersion  in  a  Sandwich  Plate  with  Low-Velocity  Interface  Layers 

Using  the  numerical  technique  described  in  the  theory  section,  we  analyzed  the  dispersion 
of  elastic  waves  in  the  five-layer  plate  shown  in  Figure  3.  The  displacement  is  assumed  to  be 
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u  (x.z.t)  =  f  (z)el'cx'*wt 


(IS) 


ll  :i  1 1  (he  layers  are  transversely  isotropic  with  the  symmetry  axes  parallel  or  perpendicular 
to  the  direction  of  wave  propagation  (x-axis),  then  the  equations  governing  the  y-componcnl  of 
the  displacement.  uv.  uncouple  from  those  governing  the  x  and  z  components.  The  former 
correspond  to  the  SI l  motion  and  the  latter  to  plane  strain  motion. 


Figures  11  and  12  show  the  disperison  curves  for  the  plane-strain  motion,  and  Figure  13 
shows  those  for  the  Sll  motion.  In  these  figures  the  vertical  axis  corresponds  10  H  = 


wh.irvMT  and  the  horizontal  axis  to  f 


kh/ir . 


Here  h  is  the  thickness  of  each  of  the  two 


low  velocity  lavers  and  C, 


is  the  shear  wave  velocity  m  these  layers.  The  properties 


of  the  stiff  layers  are  taken  to  be  =  S. 76x10’  N/m\  p' =  l  771  g/cm\  v'  =  0  3.  Those  of  the 

soft  layers  are  m  =  l.TTxlt)9  N/m!,  p  =  1.2  g/cm5.  »  =  0.3.  h/H  is  taken  to  be  0.1 

Figure  11  shows  the  first  nine  modes  for  small  values  of  tt  and  f  It  is  seen  that  these 
dispersion  curves  look  qualitatively  very  similar  to  those  for  an  isotropic  plate  However,  the 

important  difference  is  that  cut-off  frequencies  of  higher  modes  are  lowered  signfi'cantly.  Figure 
12  shows  the  first  three  modes  over  a  wide  range  of  frequency  and  wave  number.  It  is  found 
that  at  short  wavelengths  the  phase  velocity  of  the  first  two  modes  departs  significantly  from  the 
Rayleigh  wave  velocity  the  stiff  layer.  Dispersion  curves  for  SH-motion  (Figure  13)  also  show 
similar  features  The  departure  depends  on  the  ratios  of  the  elastic  properties  and  of  the 
thicknesses  of  the  low-velocity  layers  and  the  laminae.  This  effect  and  the  lowering  of  cutoff 
frequencies  should  lead  to  ultrasonic  characterization  of  interfaces. 

CONCLUSIONS 

We  have  shown  that  modeling  and  experimental  observations  of  particle-reinforced  and 
fiber-reinforced  materials  lead  to  property  characterization  of  the  reinforcing  phases.  Also,  we 
have  presented  model  calculations  of  interface  effects  on  phase  velocities  and  attenuation  of  waves 
in  a  composite  medium.  It  is  shown  that,  for  the  particular  systems  considered,  the  presence  of 
low- velocity  interface  layers  decreases  the  phase  velocities  as  well  as  the  attenuation. 
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Table  1.  Measured  and  calculated  elastic  constants  for  graphite-fiber-reinforced 
composite  and  calculated  Araphite-f iber  elastic  constants.  Except  for  dimension¬ 
less  v^,  units  are  GPa. 

22  Composite  Fiber  Composite  Composite- 
Fiber  Measured  Calculated  Calculated  Calc./Meas. 
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Figure  1.  An  oblate  spheroidal  inclusion  with  interface  layer. 

Figure  2.  A  spherical  inclusion  with  interface  layer. 

Figure  3.  A  layered  plate  with  thin  interface  layers. 


Figure  4.  For  cast-iron.  Young's  modulus  versus  graphite-particle  aspect  ratio.  Symbols 

r  represent  measurements.  Curves  represent  model-calculation  results  for  two  volume 

!  fractions:  0  10  and  0.12.  Upper,  nearly  horizontal,  curves  represent  graphite's  upper 

I  third-order-bound  (Kroner-bound)  elastic  constants.  Lower  curves  represent 

,  graphite's  lower  third-order  bound. 


Figure  5. 
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Figure  7. 

Figure  8 


Craphite-Ai  microstructure.  Transverse  section  of  parallel  7-/zm-diametcr  graphite 
fibers  distributed  nearly  homogeneously  in  aluminum  matrix.  The  fiber  volume 
fraction  equals  70  percent.  The  white  network  represents  aluminum  boundary 
regions  between  fiber  bundles  used  in  manufacture. 

For  a  composite  consisting  of  70-vol.-pct.  uniaxial  graphite  fibers  in  an  aluminum 
matrix,  the  variation  with  temperature  of  the  principal  C,j  elastic  stiffnesses.  Lines 
at  the  right  show  the  variation  of  the  aluminum-matrix  longitudinal  and  shear 
moduli,  and  G. 

Attenuation  of  a  plane  longitudinal  wave  in  a  particle-reinforced  composite  with  and 
without  interface  layers,  c  denotes  volume  function  of  inclusions. 

Attenuation  of  a  plane  shear  wave  in  a  particle-reinforced  composite  with  and 
without  interface  layers. 


Figure  9.  Phase  velocity  of  a  plane  longitudinal  wave  in  particle-reinforced  composite  with  and 
without  interface  layers. 


Figure  10.  Phase  velocity  of  a  plane  shear  wave  in  a  particle- reinforced  composite  with  and 
without  interface  layers. 


Figure  11.  Dispersion  of  Lamb  waves  in  a  sandwich  plate  at  low  frequencies. 

Figure  12.  Dispersion  of  Lamb  waves  in  a  sandwich  plate  at  finite  frequencies. 


Figure  13.  Dispersion  of  SH  waves  in  a  sandwich  plate. 
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